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1 Introduction 

Traditionally, B physics has been the domain of e^e~ machines operating on the 
T(4S') resonance or the Z° pole. But the UAl Collaboration has already shown that 
B physics is feasible at a hadron collider environment (see for example Ref. [1]). The 
first signal of fully reconstructed B mesons at a hadron collider has been published by 
the CDF Collaboration in 1992 [2]. CDF reconstructed a handful of B+ J li)K^ 
events in a data sample of 2.6 pb~^ taken during the Tevatron Run at the end of the 
1980's. Since then experimental techniques improved significantly. Especially with 
the development of high precision silicon vertex detectors, the study of B hadrons is 
now an established part of the physics program at hadron colliders. 

The CDF and DO experiments can look back to an already successful B physics 
program during the 1992-1996 Run I data taking period (for a review of B physics 
results from, for example, CDF in Run I see Ref. [3]). Nowadays, B physics results 
from a hadron collider are fully competitive with the e^e~ B factories. As discussed 
later in this review, with the operation of a hadronic track trigger, CDF reconstructs 
fully hadronic B decay modes without leptons in the final state. In many cases, the 
measurements performed at the Tevatron Collider are complementary to the B fac- 
tories. For example, no 5° mesons or baryons containing h quarks are produced on 
the T(4S') resonance. 

B hadrons not produced at the B factories are the topics of this review. We discuss 
the spectroscopy of excited B states (5**, S**) and the observation of the Sf, baryon 
at the Tevatron. The second part of this review discusses the decays of B hadrons 
and measurements of branching fractions. We focus on charmless two-body decays of 
B — > h'^h~. We end this article by summarizing our finding in the conclusions. 
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Figure 1: Tevatron (a) initial store luminosity from 2002-2006 and (b) delivered 
luminosity per calendar year. 



2 The Tevatron with the CDF & DO Experiments 

The Fermilab accelerator complex has undergone a major upgrade in preparation 
for Tevatron Run II. The centre-of-mass energy has been increased to 1.96 TeV as 
compared to 1.8 TeV during Run I and the Main Injector, a new 150 GeV proton 
storage ring, has replaced the Main Ring as injector of protons and anti-protons into 
the Tevatron. The present bunch crossing time is 396 ns with a 36 x 36 pp bunch 
operation. The luminous region of the Tevatron beam has an RMS of ~ 30 cm along 
the beamline (2;-direction) with a transverse beamwidth of about 25-30 /im. 

The initial Tevatron luminosity steadily increased from 2002 to 2006 as shown in 
Figure 1(a). By the end of 2006, the peak luminosity reached by the Tevatron is 
> 25 ■ 10'^^ cm~^s~^. The increase in accelerator performance throughout Run II can 
also be seen by the delivered luminosity per calendar year as displayed in Figure 1(b). 
The total integrated luminosity delivered by the Tevatron to CDF and DO at the end 
of 2006 is ~ 2.2 fb~^ with about 1.8 fb~^ recorded to tape by each the CDF and 
DO experiments. However, most results presented in this review use about 1 fb^^ of 
data. 

The CDF detector improvements for Run II [4] were motivated by the shorter ac- 
celerator bunch spacing and the increase in luminosity by an order of magnitude. All 
front-end and trigger electronics has been significantly redesigned and replaced. A 
DAQ upgrade allows the operation of a pipelined trigger system. CDF's tracking sys- 
tem was completely renewed for Run II. It consists of a Central Outer Tracker (COT) 
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with 30200 sense wires arranged in 96 layers, between 40 and 137 cm in radius, 
organized into eight alternating axial and ±2° stereo super-layers. The transverse 
momentum resolution is o-p^/pr — 0.15%pr/(GcV/c). The specific energy loss by 
ionization (dE/dx) of charged particles in the COT is measured from the amount of 
charge collected by each wire. The Run II silicon vertex detector consists of seven 
double sided layers and one single sided layer mounted on the beam pipe covering a 
total radial area from 1.5-28 cm. The silicon vertex detector covers the full Tevatron 
luminous region and allows for standalone silicon tracking up to a pseudo-rapidity \r]\ 
of 2. The forward calorimeters have been replaced by a new scintillator tile based 
plug calorimeter which gives good electron identification up to |?7| = 2. The upgrades 
to the muon system almost double the central muon coverage and extend it up to 
\r]\ ~ 1.5. The most important improvements for B physics in Run II are a Silicon 
Vertex Trigger (SVT) and a Time-of-Flight (ToF) system with a resolution of about 
100 ps. The later employs 216 three-meter-long scintillator bars located between the 
outer radius of the COT and the superconducting solenoid. The Time-of-Flight sys- 
tem is most beneficiary for the identification of kaons with a 2 cr-separation between 
TT and K ior p < 1.6 GeV/c. 

The DO detector also went through a major upgrade before the beginning of 
Run II [5]. The inner tracking system was completely replaced and includes a new 
Silicon tracker surrounded by a Scintillating Fiber tracker, both of which are en- 
closed in a 2 Tesla solenoidal magnetic field. Pre-shower counters are located before 
the uranium/liquid-argon calorimeter to improve the electron and photon identifica- 
tion. The already excellent muon system has been further improved by adding more 
shielding to reduce beam background. The Run II DO detector has excellent tracking 
and lepton acceptance. Tracks with pseudo-rapidity as large as 2.5-3.0 {9 ~ 10°) 
and transverse momentum px as low as 180 MeV/c can be reconstructed. The muon 
system can identify muons within I77I < 2.0. The minimum pr of the reconstructed 
muons varies as a function of 77. In most of the results presented, muons were required 
to have pr > 2 GeV/c. 

2.1 Triggers for B Physics 

The total inelastic pp cross section at the Tevatron is about three orders of magnitude 
larger than the b quark production cross section. The CDF and DO trigger system 
is therefore the most important tool for finding B decay products. In addition, the 
cross section for b quark production is steeply falling. It drops by almost two orders 
of magnitude between a b quark pt of about 8 GeV/c and 25 GeV/c. To find B decay 
products in hadronic collisions, it is desirable to go as low as possible in the decay 
products transverse momentum, exploiting as much as possible of the steeply falling 
b cross section. Of course, the limiting factor is the bandwidth of the experiment's 
data acquisition system. 
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In Run I, all B physics triggers at CDF and DO were based on leptons including 
single and dilepton triggers. In Run II, both experiments still exploit heavy flavour 
decays which have leptons in the final state. Identification of dimuon events down 
to very low momentum is possible, allowing for efficient J/ip ^ ^'^fi~ triggers. As a 
consequence, both experiments are able to fully reconstruct B decay modes involving 
J/ip^s. Both experiments also use inclusive lepton triggers designed to accept semilep- 
tonic B ii/^X decays. DO has an inclusive muon trigger with excellent acceptance, 
allowing them to accumulate very large samples of semileptonic decays. The CDF 
semileptonic triggers require an additional displaced track associated with the lepton, 
providing cleaner samples with smaller yields. 

In addition, the CDF detector has the ability to select events based upon track 
impact parameter. The Silicon Vertex Trigger gives CDF access to purely hadronic 
B decays and makes CDF's B physics program fully competitive with the one at 
the e~^e~ B factories. The hadronic track trigger is the first of its kind operating 
successfully at a hadron collider. It works as follows: With a fast track trigger 
at Levell, CDF finds track pairs in the COT with pr > 1.5 GeV/c. At Level 2, 
these tracks are linked into the silicon vertex detector and cuts on the track impact 
parameter (e.g. d > 100 //m) are applied. The SVT track impact parameter resolution 
is about 50 fim including a 33 /im contribution from the transverse beam spreading. 
The original motivation for CDF's hadronic track trigger was to select B^ — > tttt 
decays to be used for CP violation studies. With the different B trigger strategies 
above, the Collider experiments are able to trigger and reconstruct large samples of 
heavy flavour hadrons. 

3 Spectroscopy 

3.1 Study of Orbit ally Excited B Mesons 

The spectroscopy of excited meson states containing b quarks is not well studied. 
Only the stable 0~ ground states B~^, B^ and B^ and the excited 1~ state B* are 
established [6]. Quark models predict the existence of two wide {Bq and Bl) and two 
narrow (5° and i?2*) bound P-states [7]. The wide states decay through an 5*- wave 
and therefore have a large width of a couple of hundred MeV/c^, which makes it 
difficult to distinguish such states from combinatoric background. The narrow states 
decay through a D-wave (L = 2) and thus should have a small width of around 
1 MeV/c^ [8, 9]. Almost all previous observations [10, 11] of the narrow P-states 
Bl and B2* have been made indirectly using inclusive or semi-exclusive B decays 
which prevented the separation of both states and a precise measurement of their 
properties. In contrast, the masses, widths and decay branching fractions of these 
states are predicted with good precision by the theoretical models [8, 9]. 
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Figure 2: Result of the fit to the B** mass difference (a) Am = m{B7c) — m{B) from 
DO and (b) Q = m(i?7r) — m{B) — m(7r) from CDF in the B^ J /ipK^ channel 
and (c) in the B^ D^tc'^ mode. 



B^ and B2* candidates are reconstructed in the following decay modes: B^ 
B*+TT- with B*+ and B^* B*+tt- with B*+ 5+7 as well as ^ 

B~^7T^. In both cases the soft photon from the B* decay is not reconstructed result- 
ing m a shift of about 46 MeV/c^ in the mass spectrum. DO reconstructs the B^ 
candidates in the fully reconstructed mode B^ J/ipK^ with Jjip ^ while 
CDF selects B^ mesons in addition through the B^ D^n^ mode with K^tc^. 
The CDF analysis is based on 360 pb~^ of data resulting in a B^ —>■ J/tpK^ signal 
of 1867 ± 64 events and 2182 ± 54 candidates in the B~^ — * D^n^ channel. The DO 
measurement employs 1 fb~^ of Run II data and finds a signal peak of 16 219 ± 180 
events attributed to the decay B^ J/ipK'^. 

DO presents their measured mass distribution as Am = m{BiT) — m{B) as shown 
in Figure 2(a), while CDF plots Q = m{B7r) — m{B) — m(7r) as displayed in Fig. 2(b) 
and (c). Clear signals for the narrow excited B states are observed: CDF reconstructs 
80 ±18 events in B^ J/ipK^ and 106 ±20 events in the B^ D^7t~^ channel while 
DO observes a total of 504 ± 80 candidates for the narrow B** states. The measured 
masses are reported as m{B\) = 5720.8 ± 2.5 ± 5.3 MeV/c^ and m{B2*) — 'm{B\) = 
25.2 ± 3.0 ± 1.1 MeV/c^ from DO, while CDF quotes m(5°) = 5734 ± 3 ± 2 MeV/c^ 
and m(i?2*) = 5738±5± 1 MeV/c^. Clearly these preliminary results are not in good 
agreement. CDF currently works on an update of their analysis using 1 fb~^ of data. 
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3.2 Observation of Orbit ally Excited Bsj Mesons 

The properties of (bs) excited meson states and the comparison with properties of 
excited states in the {Jm) and {bd) systems provide good tests of various models 
of quark bound states. These models [7, 8, 12] predict the existence of two wide 
resonances {B*q and B*^) and two narrow and ) bound P-states. The wide 
states decay through an 5'-wave and therefore have a large width of a couple of 
hundred MeV/c^. This makes it difficult to distinguish such states from combinatoric 
background. The narrow states decay through a D-wave (L = 2) and therefore should 
have a small width of around 10 MeV/c^ [9]. If the mass of the Bgj (J = 1,2) is 
large enough, then the main decay channel should be B^*^K as the B^ti decay mode 
is not allowed by isospin conservation. Previous observations [10] of the narrow Bgj 
P-states have been made indirectly preventing the separation of both states. 

Bg^ and candidates are reconstructed in the following decay modes: B^^ — > 
B*+K- with B*+ B+-f and B^* B*+K- with B*+ B+-f as well as B°; 
B^K~ . In both cases the soft photon from the B* decay is not reconstructed resulting 
in a shift in the mass spectrum. DO reconstructs the P+ candidates in the fully 
reconstructed mode B^ —>■ J/ipK^ with J/tp ^ l^^ while CDF selects B^ mesons 
in addition through the P+ D\+ mode with P)° K'-k^. The CDF and DO 
measurements are each based on 1 fb~^ of Run II data. The CDF analysis finds 
~ 31 000 S+ ^ J/ipK"^ events and ~ 27200 candidates in the P+ ^ P'°7r+ channel. 
DO uses a signal of 16 219 ± 180 B^ events from the decay B^ J/ipK'^. Both 
experiments present their measured mass distribution in the quantity Q = m{BK) — 
m{B) — m{K) as displayed in Figure 3(a) and (b). 

A clear signal at Q ~ 67 MeV/c^ is observed by CDF and DO (see Fig. 3), which 
is interpreted as the state. CDF reconstructs 95 ± 23 events in the peak at 
Q = 67.0 MeV/c2 while DO reports 135 ± 31 events at Q = 66.4 ± 1.4 MeV/c^. 
In addition, CDF observes 36 ± 9 events in a peak at Q ~ 10.7 MeV/c^ which is 
interpreted as first evidence for the B^^ state. The measured masses are reported 
as m(P°2*) = 5839.1 ± 1.4 ± 1.5 MeV/c^ from DO, while CDF quotes m(P°i) = 
5829.4 ± 0.2 ± 0.6 MeV/c^ and m(P°2*) = 5839.6 ± 0.4 ± 0.5 MeV/c^ The results 
from CDF and DO are in good agreement. 

3.3 Observation of E;, Baryons 

Until recently only one bottom baryon, the A°, has been directly observed. At present 
the CDF collaboration has accumulated the world's largest data sample of bottom 
baryons, due to a combination of two factors - the CDF displaced track trigger, and 
the ~ 1 fb~^ of integrated luminosity delivered by the Tevatron. Using a sample of 
fully reconstructed A° — > A+tt" candidates collected with the displaced track trigger, 
CDF searched for the decay eJ*^^ A^tt^. 
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Figure 3: Result of the fit to the B*j mass difference Q = m{BK) — m{B) — m{K) 
from (a) CDF and (b) DO. 



The QCD treatment of quark-quark interactions significantly simplifies if one of 
the participating quarks is much heavier than the QCD confinement scale Aqcd- 
In the limit of uiq — > oo, where mq is the mass of the heavy quark, the angular 
momentum and flavour of the light quark become good quantum numbers. This 
approach, known as Heavy Quark Effective Theory (HQET), thus views a baryon 
made out of one heavy quark and two light quarks as consisting of a heavy static 
color field surrounded by a cloud corresponding to the light diquark system. In 
SU(3) the two quarks are in diquark form 3 and 6 according to the decomposition 
3 ® 3 = 3 © 6, leading to a generic scheme of baryon classification. Diquark states 
containing quarks in an antisymmetric flavour configuration, [gi, 52], are called A- type 
whereas states with diquarks containing quarks in a flavour symmetric state, {gi, ^2}, 
are called S-type. 

In the S-type ground state the light diquark system has isospin 1 = 1 and = 1^. 
Together with the heavy quark this leads to a doublet of baryons with = (Sf,) 

and = I"*" (S^). The ground state E-type baryons decay strongly to A-type 
baryons by emitting pions. In the limit ttiq 00, the spin doublet {Sft, S^} would 
be exactly degenerate since an infinitely heavy quark does not have a spin interaction 
with a light diquark system. As the heavy quark is not infinitely massive, there 
will be a small mass splitting between the doublet states and there is an additional 
isospin splitting between the and states [13]. There exist a number of 

predictions for the masses and isospin splittings of these states using HQET, non- 
relativistic and relativistic potential models, l/N^ expansion, sum rules and lattice 
QCD. References [13, 14] contain some of the existing theoretical estimates, while 
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Eft property 


Expected value [MeV/c^] 


m(Eft) - m(A^) 
m(E,*) - m(E5) 
m(Eft--) - m(E+) 
r(Sft), r(E^) 


180 - 210 

10 - 40 
5 - 7 
-8, -15 



Table 1: General range of theoretical predictions for the E^*^^ states from Refer- 
ences [13, 14]. 



Table 1 summarizes the range of predictions. The natural width of E^ baryons is 
expected to be dominated by single pion transitions. Decays of the type I^c,b Ac,b1 
are expected to have significantly smaller (~ 100 keV/c^) partial widths than the 
single pion transition, and are thus negligible. The partial width of the P-wave one- 
pion transition thus depends on the available phase space. 

In analogy with the B meson hadronization chain, in this analysis events are 
separated into "same charge" or SC and "opposite charge" or OC combinations. As 
the is neutral, the charge of the soft pion track determines the charge of the 
E5 baryon, and there will be E^ signals for both positive and negative pions. SC 
(OC) is defined as events where the E5 pion has the same (opposite) charge as the 
pion from the A5 decay. With these definitions, the SC distribution contains all E^*'' 

and E^*"* candidates while OC contains E^*-*^ and E^*^"*". 

The present analysis is based on events collected by the CDF detector from 2002 
through February 2006, with an integrated luminosity oi C — 1070 ±60 pb~^. Events 

collected on the two track trigger are used to reconstruct the decay chain A° — >■ 
A^TT, A^ pK^n^. CDF reconstructs a A^ yield of approximately 2800 candidates 
in the signal region m(A^) e [5.565,5.670] GeV/c^, with the A^ mass plot shown in 
Figure 4. 

To separate out the resolution on the mass of each A° candidate, CDF searches for 
narrow resonances in the mass difference distribTition of Q = m(A°7r) — m{A^) — m.^. 
Unless explicitly stated, E5 refers to both the J = ^ (^6^) ^^'^ ~ I (^ti^) states. 
There is no transverse momentum cut applied to the pion from the E^ decay, since 
these tracks are expected to be very soft. In order to perform an unbiased search, the 
cuts for the E5 reconstruction are optimized first with the E^ signal region blinded. 
From theoretical predictions the E5 signal region is chosen as 30 < Q < 100 MeV/c^, 
while the upper and lower sideband regions of < Q < 30 MeV/c^ and 100 < 
Q < 500 MeV/c^ represent the E^ background. The signal for the optimization is 
taken from a PYTHIA Monte Carlo E;, sample, with the decays E^ — > A°7r, A° 
A+7r~, A+ — ^pX^TT"^ forced. 

The backgrounds under the A° signal region in the A° mass distribution will 
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Figure 4: Fit to the invariant mass of — * A+vr candidates. The sohd blue hne is 
the total fit, while the primary background sources are listed in the legend. 



also be present in the Q-distribution. The primary sources of background are 
A° hadronization and underlying event, hadronization and underlying event of other 
B meson reflections and combinatorial background underneath the A° peak. The 
percentage of each background component in the A^ signal region is derived from the 
A° mass fit, and is determined as 86% A° signal, 9% backgrounds and 5% combina- 
torial background. Other backgrounds {e.g. from 5-track decays where one track is 
taken as the tts^ candidate) are negligible, as confirmed in inclusive single-6-hadron 
Monte Carlo samples. 

Upon unblinding the Q signal region, there is an excess observed in data over 
predicted backgrounds. The excess over background is shown in Table 2. CDF 
performs a simultaneous unbinned likelihood fit to SC and OC data. To the already 
described background components, four peaks are added, one for each of the expected 
Eft states. Each peak is a sum of two Breit-Wigner shapes, each convoluted with two 
Gaussian resolution functions. The detector resolution has a dominant narrow core 
and a small broader shape describing the tails where the PDF for each peak takes 
both into account. Due to low statistics, CDF constrains m(S^^) — m(E^) and 
m(E^~) — m(S^) to be the same. The results of the fit are given in Tab. 3 and 
displayed in Fig. 5(a). 

All systematic uncertainties on the mass difference measurements are small com- 
pared to their statistical errors. The systematic errors from the tracking sources are 
determined by comparing the mean and the width of the peak in m{D*^) — m{D'^) 
between data and Monte Carlo simulation split up in several regions of track px- The 
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Sample 


Data events 


Bkg events 


Data excess over bkg 


Same charge 
Opposite charge 


416 
406 


268 
298 


148 
108 



Table 2: Summary of the number of events in the Q signal region {Q e [0.03, 
0.1] GeV/c^) for data and predicted background. 



Parameter 


Value 


Parabolic Error 


MINOS Errors 


Q(S+) (MeV/c^) 
g(S,-) (McV/c2) 
g(S*) - Q(Efe) (MeV/c2) 


48.4 
55.9 
21.3 


2.02 
0.963 
1.93 


(+2.02, -2.29) 
(+0.990, -0.959) 
(+2.03, -1. 94) 


events 


29 


12.0 


(+12.4, -11.6) 


events 


60 


14.3 


(+14.8, -13.8) 


E^^ events 


74 


16.8 


(+17.2, -16.3) 


E^^ events 


74 


17.8 


(+18.2, -17.4) 


-ln(Likelihood) 


-24553.5 







Table 3: Fit parameters and error values from the fit to data. Positive and negative 
errors are quoted separately as the error range is asymmetric. 



largest discrepancy of the D*^ peak is 0.06 MeV/c^ which is taken as the systematic 
error for all four peaks. The discrepancy in the mass resolution could be as large 
as 20%. The effect of a broader resolution is evaluated via a sample of Toy Monte 
Carlo experiments. The remaining systcmatics come from assumptions made in the 
fit to the data, such as the use of fixed background shapes. For the parameters as- 
sociated with an individual systematic uncertainty. Toy MC samples are generated 
where these parameters are varied. The sample is then fit with both the default fit 
and the fit with varied parameters. The difference between fit parameter values in 
the varied fit and the default fit is caused by the systematic variation and constitutes 
the associated systematic error. 

To evaluate the significance of the measurement, the null hypothesis is tested. 
The data is fit with no signal and with the standard fit using four peaks. Then 
the likelihood ratio is computed as LR = L1/L2, where L2 is the four signal peak 
hypothesis and Li is the corresponding hypothesis with no peaks. The result of this 
fit is shown in Figure 5(b) and a likelihood ratio of ~ 10~^^ is obtained indicating the 
observation of the E^*-*^ states. 

To summarize, the lowest lying charged A°7r resonant states are observed in 1 fb~^ 
of data collected by the CDF detector. These are consistent with the lowest lying 
charged E^*^^ baryons. The Q values of E^ and E^, and the E^-E^ mass difference. 
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Q = m(A°7r) - m(A°) - (GeV/c^) Q = m(A°7r) - m(A°) - (GeV/c^) 



Figure 5: (a) Simultaneous fit to the states and (b) with an alternate signal 
description assuming no signal is present (null hypothesis). 



are measured to be: 

- m(Sfe-) - m(A°) - m(7r) = 55.9 ± 1.0 (stat) ±0.1 (syst) MeV/c^, 

- m(E^) - m(Ag) - m(7r) = 48.4t|^ (stat) ±0.1 (syst) MeV/c^, 

- m(Er) - m(Sfc-) = m(S^+) - m(S^) = 21.3;?;^ (stat) (syst) MeV/c^. 

Using the best CDF mass measurement for the mass, which is m{A^) = 5619. 7± 
1.2 (stat) ±1.2 (syst) MeV/c^, the absolute mass values and number of events are: 

- m(S^) = 5808tli (stat) ±1.7 (syst) MeV/c^, N(S^) = 29tni (stat) (syst), 

- m(Eft-) = 5816tl:[j (stat) ±1.7 (syst) MeV/c^, N(S^) = 60tl|^ (stat) (syst), 

- m(E*+) = 5829tJ| (stat) ±1.7 (syst) MeV/c^, N(S*+) = 74tJ^-2 (stat) ±1^-^^ (syst), 

- m(E*-) = 5837;?-^ (stat) ±1.7 (syst) MeV/c^, N(S*-) = 74tJ^i (stat) I^^q^ (syst). 



4 Decay of S Hadrons 

In this Section we focus on a new CDF result involving the branching fractions and 
time-integrated direct CP asymmetries for and decay modes into pairs of 
charmless charged hadrons B h^h~ . 
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4.1 Results from Charmless Two-Body Decays B h^h 

The decay modes of B mesons into pairs of charmless pseudo-scalar mesons arc effec- 
tive probes of the quark-mixing matrix (CKM) and sensitive to potential new physics 
effects. The large production rate of B hadrons at the Tevatron allows measuring 
such decays in new modes, which arc important to supplement our understanding of 
B meson decays. The still unobserved 5^ K^n^ decay mode could be used to mea- 
sure the angle 7 [15] of the CKM unitarity triangle and its CP asymmetry could be a 
powerful model- independent test of the source of direct CP violation in the B meson 
system [16]. This may provide useful information to solve the current discrepancy 
between the asymmetries observed in the neutral and charged B modes [17]. The 
B^ 7r"'"7r~ and B^ K^K^ decay channels proceed only through annihilation 
diagrams, which are currently poorly known and constitute a source of significant 
uncertainty in many theoretical calculations [18, 19]. A measurement of both modes 
would allow a determination of the strength of penguin-annihilation diagrams [20] . 

4.1.1 Data Selection 

CDF analysed a sample (integrated luminosity £ ~ 1 fb~^) of pairs of oppositely 
charged particles with pr > 2 GeV/c and Pt(1) +Pt(2) > 5.5 GeV/c, used to form 
B^^^ meson candidates. In addition, the trigger required a transverse opening-angle 
20° < A(j) < 135° between the two tracks, to reject background from particle pairs 
within the same jet and from back-to-back jets. In addition, both charged particles 
are required to originate from a displaced vertex with a large impact parameter do 
(100 //m < do < 1 mm), while the B^^^ meson candidate is required to be produced 
in the primary pp interaction {do{B) < 140 /im) and to have traveled a transverse 
distance L^y^B) > 200 /um. 

In the offline analysis, an unbiased optimization procedure determines a tightened 
selection on track-pairs fit to a common decay-vertex. CDF chooses selection cuts 
minimizing directly the expected uncertainty (through several pseudo-experiments) of 
the physics observables to be measured. CDF decided to use two different sets of cuts, 
optimizing separately the measurements of Acp{B^ — > K^n') and i3(-B° — > K~t{^). 
For the latter, the sensitivity for discovery and limit setting [21] was optimized rather 
than the statistical uncertainty on the particular observational parameter, since this 
mode had not yet been observed. It is verified that the former set of cuts is also 
adequate to measure other decay rates of the larger yield modes {B^ tt'^tt^ , B^ 
K~^K~), while the latter, tighter set of cuts, is well suited to measure the decay 
rates and CP asymmetries related to the rare modes {B^ tt+tt", B° — > K'^K~, 

In addition to tightening the trigger cuts in the offline analysis, other discrim- 
inating variables such as the isolation of the B^^^ meson and the information pro- 
vided by the 3D reconstruction capability of the CDF tracking system are used. 
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Figure 6: (a) Invariant mass distribution of B ^ h^h^ candidates passing all selection 
requirements optimized to measure ^3(5° K^tt^), using the pion mass assumption 
for both decay products. The cumulative projections of the likelihood fit for each 
mode are overlaid in (b). 



allowing a great improvement in the signal purity. Isolation is defined as I{B) = 
Pt{B)/[pt{B) + J2iPT{i)], in which the sum runs over every other track within a cone 
of radius one in the rj — (p space around the 5°^^ meson flight-direction. By requiring 
I{B) > 0.5 the background is reduced by a factor four while keeping almost 80% of 
the B signal. The 3D silicon tracking allows to resolve multiple vertices along the 
beam direction and to reject fake tracks reducing the background by another factor 
of two, with small inefficiency on the signal. The resulting vrvr invariant mass distri- 
bution shown in Figure 6(a) display a clean signal of B ^ h^h^ decays. In spite of a 
good mass resolution 22 MeV/c^), the various B h^h^ modes overlap into an 
unresolved mass peak. 

4.1.2 Fit of Sample Composition 

The resolution in invariant mass and in particle identification is not sufficient for sep- 
arating individual decay modes on an event-by-event basis. Therefore CDF performs 
an unbinned maximum likelihood fit, combining kinematic and particle identifica- 
tion information, to statistically determine the contribution of each mode and the 
CP asymmetries. For the kinematic portion, CDF uses three loosely correlated ob- 
servables to summarize the information carried by all possible values of invariant mass 
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of the B candidate, resulting in different mass assignments to the two outgoing parti- 
cles. These are: (a) the mass m^^ calculated with the charged pion mass assignment 
to both particles, (b) the signed momentum imbalance a = (1 — Pi/p2)qi, where pi 
{P2) is the lower (higher) of the particle momenta, and qi is the sign of the charge 
of the particle of momentum pi, and (c) the scalar sum of the particle momenta 
Ptot = Pi + P2- Using these three variables, the mass of any particular mode 71112 can 
be written as: 



= rnl^ — 2ml + rn\ + m\- 2Jpl + mlJpl + ml - 2Jpl + mjJpl + ml, (1 



Pi = 7^ r^Ptot, P2 = - ^Ptou (2) 

2 — |q;| 2 — \a\ 

where mi (m2) is the mass of the lower (higher) momentum particle. For simplicity 
Eq. (1) is written as a function of pi and p2 instead of a and ptot but in the likelihood 
fit it is used as a function of a and ptot- 

Particle identification (PID) information is summarized by a single observable k 
for each track defined as 

_ dE/dx - dE/dx{7r) 
~ dE/dx{K) - dE/dx{7r) ■ 

With the chosen observables, the likelihood contribution of the i*'* event is written 
as: 

A = (1 - &) E f,^^"^ + b {fA^Tj^T + (1 - /a)>C|^-£^^°) (4) 



where: 



= R{m^n - Mj{a,ptot), a, Ptot) Pj {a, Ptot), (5) 
>Ca™ = A(m^^|c2, mo)PA(Q;,Ptot), (6) 
4'" = e^^"^"''^E(«,Ptot), (7) 

^m^) = ^i(E,A)(Ki, K2, a, Ptot)- (8) 

The index 'A(E)' labels the physical (combinatorial) background-related quantities, 
the index j runs over the twelve distinguishable B — > h'^h~ and A° — > ph modes 
(Fig. 7), and fj are their respective fractions, to be determined by the fit together with 
the total background fraction b and with the fraction of the physical (combinatorial) 
background fA(E)- The conditional probability density Rlm^^j^ — Aij{a,ptot),(^,Ptot) 
is the mass resolution function of each mode j when the correct mass is assigned to 
both tracks. In fact, the average mass M.j{a,ptot) is the value of m^^^^ obtained from 
Eq. (1) by setting the appropriate particle masses for each decay mode j. Making a 
simple variable change, R^mj^T, — Aij{a,ptot),a,Ptot) = Rijrij — m^o(^QQ^j^o-^,a,ptot) is 
obtained where mj is the invariant mass computed with the correct mass assignment 
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Figure 7: Average mTTTr vs a for simulated samples of (a) and (b) A^ candidates, 
where self-tagging final states {K'^7r~ and K~Tr'^, ph~ and ph'^) are treated sepa- 
rately. The corresponding plots for the are similar to B^ but shifted for the mass 
difference. 



to both particles for each mode j. R is parameterized using the detailed detector 
simulation [22]. To take into account non-Gaussian tails due to the emission of 
photons in the final state, CDF includes in the simulation soft photon emission of 
particles in agreement with recent QED calculations [23]. CDF checks the quality 
of the mass resolution model using about 500K D^^ — > K^n^ decays as shown in 
Figure 8(a). The mass hne-shape of the — > K~n'^ peak is fitted fixing the signal 
shape from the model, only allowing to vary the background function. CDF obtains 
good agreement between data and simulation. In Eq. (5) the nominal B^, B^ and 

masses as measured by CDF [24] are used in order to cancel common systematic 
uncertainties. The background mass distribution is determined in the fit by varying 
the parameters ci, C2 and mo in Eq. (6,7). The probability Pj{a,ptot) is the joint 
probability distribution of {a,ptot) and is parameterized for each mode j by a product 
of polynomial and exponential functions fitted to Monte Carlo samples produced by 
a detailed detector simulation [22]. The background function Pa{e) is obtained from 
the mass sidebands of the data. 

A sample of 1.5M D^-k^ [ir~7r+]7r+ decays, where the decay prod- 

ucts are identified by the charge of the D*^ pion, was used to calibrate the dE/dx 
response over time and over the entire tracking volume, and to determine the F func- 
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Figure 8: Tagged D° K-7i+ decays from D*+ D°7r+ [K-n+]n+. (a) Check of 
the mass hne shape template performing a 1-dimensional binned fit where the signal 
mass line shape is completely fixed from the model, (b) Distribution of dE/dx (mean 
COT pulse-width) around the average pion response for calibration samples of kaons 
(left) and pions (right). 



tions in Eq. (8). Using a > 95% pure D'^ sample, CDF obtains a 1.4 a separation 
between kaons and pions as shown in Fig. 8(b), corresponding to an uncertainty on 
the measured fraction of each class of particles that is just 60% worse than the uncer- 
tainty attainable with ideal separation. The background term in Eq. (8) is similar to 
the signal terms, but allows for independent pion, kaon, proton, and electron compo- 
nents, which are free to vary independently. Muons are indistinguishable from pions 
with the available dE/dx resolution. 

4.1.3 Fit Results 

CDF performs two separate fits. The first one uses the cuts optimized to measure 
the direct Acp{B^ i^+vr") and the second one is optimized to measure B^B'^ — > 
K~7r~^). Significant signals are seen for the —>■ tt+tt", 5° — >■ K~^7r~, and B^ —>■ 
K^K~ modes, previously observed by CDF [25]. Three new rare modes are observed 
for the first time: 5° K~Tr~^, pn^ and — > pK~ , while no evidence is 

obtained for the B^ —>■ tt+tt^ and B^ —>■ K^K~ decay channels. 

To convert the yields returned from the fit into relative branching fractions, CDF 
applies corrections for efficiencies of trigger and offline selection requirements for the 
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different decay modes. The relative efficiency corrections between various modes 
do not exceed 20%. Most corrections are determined from the detailed detector 
simulation [22] , with some exceptions which are measured using data. A momentum- 
averaged relative isolation efficiency between and mesons of 1.07 ±0.11 is 
determined from fully-reconstructed samples of B^ — > J/ip(f), B^ —>■ D~t:~^, B^ — > 
J/ipK*'^, and B^ D~7r'^. The lower specific ionization of kaons with respect to 
pions in the drift chamber is responsible for a ~ 5% lower efficiency to reconstruct a 
kaon. This effect is measured in a sample of K~7r~^7r~^ decays triggered with the 

two track trigger, using the unbiased third track. The only correction needed by the 
direct CP asymmetries Acp{B° ^ K+tt') and Acp{B° K-7r+) is a < 0.6% shift 
due to the different probability for K'^ and K~ to interact with the tracker material. 
This correction uses a sample of IM prompt — > K~7r'^ decays reconstructed and 
selected with the same criteria as the B — > h'^h~ decays. Assuming the Standard 
Model expectation AcpiD^ K~tt^) = 0, the difference between the number of 
reconstructed D*^ K^n^ decays and D — > K~^tt~ provides a measurement of 
the detector-induced asymmetry between i^+7r~ and K^n'^ final states. Since CDF 
uses the same fit technique developed for the B — > h'^h" decays, this measurement 
provides also a robust check on all possible charge asymmetry biases of the detector 
and d£'/da; parameterizations. 

The 5° K+R- and 5° ^ tt+tt' modes require a special treatment, since they 
contain a superposition of the fiavour eigenstates of the B^ meson. Their time evolu- 
tion might differ from the one of the flavour-specific modes if the width difference AF^ 
between the S° mass eigenstates is signiflcant. The current result is derived under 
the assumption that both modes are dominated by the short-lived B^ component, 
that means F^ = F^^, and AF^/F^ = 0.12 ± 0.06 [26, 27]. The latter uncertainty is 
included in estimating the overall systematic uncertainty. 

The dominant contributions to the systematic uncertainty are as follows. The 
statistical uncertainty on the isolation efficiency modes), the uncertainty on the 
dE/dx calibration and parameterization and the uncertainty of the combinatorial 
background model. The first one is the larger systcmatics of all measurements with 
the meson B^ in the initial state except for Acp{B'^ — > K~tt~^). This uncertainty 
is preliminary and conservative, a significant improvement is expected for the final 
results. The second one, due to dE/dx, is a large systematics of all measurements, 
although the parameterization of the specific ionization dE/dx is very accurate. The 
fit of the sample composition is very sensitive to the FID information. The third 
systematic error is due to the statistical uncertainty of the possible combinatorial 
background models and it is a dominant systematics for the observables of the rare 
modes. Smaller systematic uncertainties are assigned for the trigger efficiencies, phys- 
ical background shapes and kinematics, and the B meson masses and lifetimes. 

The measured relative branching fractions are listed in Table 4, where fd and fg 
indicate the respective production fractions of B^ and B^ mesons from the fragmen- 
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Mode 



signal 



Quantity 



Measurement 



B [IQ- 



K+TT- 

^ TT+TT- 

B2 ^ K+K- 



4045 ± 84 
1121 ±63 
1307 ±64 



Acp{B^) 

B(B°^7i-+7r^) 
B(BO^K+-K-) 
/. BjEl^K+K-) 
fi B(B"^K+7r-) 



-0.086 ±0.023 ±0.009 
0.259 ±0.017 ±0.016 
0.324 ±0.019 ±0.041 



5.10 ±0.33 ±0.36 
24.4 ± 1.4 ±4.6 



/; B(5^^-ft— 7r+) 
h B{BO^K+-K-) 

Acp{B^) 

At{BI) 

f, BiB°^^+^-) 
h B{B"^K+w-) 



B 







K-TT 



+ 



230 ± 34 ± 16 



B2 TT+TT- 26 ± 16 ± 14 



B° K+K- 61 ± 25 ± 35 



156 ± 20 ± 11 
110 ±18 ±16 



B{B"- 



■K+K-) 



g(Ag^p7r-) 
B(h.l^pK-) 



0.066 ±0.010 ±0.010 
0.39 ±0.15 ±0.08 
-3.21 ±1.60 ±0.39 

0.007 ±0.004 ±0.005 
0.020 ± 0.008 ± 0.006 
0.66 ±0.14 ±0.08 



5.0 ± 0.75 ± 1.0 



0.53 ±0.31 ±0.40 
(< 1.36 @ 90% CL) 
0.39 ±0.16 ±0.12 
(< 0.7 @ 90% CL) 



Table 4: Results on data sample optimized to measure Acp{B^ — > i^+7r ) (top) 
and B{B^g K~7r^) (bottom). Absolute branching fractions are normalized to the 
world-average values B(S° ^ K+n') = (19.7 ± 0.6) x 10"^ /, = (10.4 ± 1.4)% and 

)-nB^-^-^^-) The first quoted 



/, = (39.8 ± 1.0)% [17]. We use Acp{B'') = ||r 



B{B''^^K+n-)+B{BO^K-K+) " ' ^ r{B^^^K+n-)-r{BO^K-K+} 

uncertainty is always statistical, the second is systematic. 



tation of b quarks in pp collisions. An upper limit is also quoted for modes in which 
no significant signal is observed [28]. The absolute branching fraction results listed 
are obtained by normalizing the data to the world-average of B{B^ K^7i~) [17]. 

CDF reports the first observation of three new rare charmless decays B^ K'n'^, 
A^ — > pn' and A° — > pK" with a significance respectively of 8.2a-, 6a and ll.Scr. 
The significance includes both statistical and systematic uncertainty. The statistical 
uncertainty to evaluate the significance is estimated using several pseudo-experiments 
with no contributions from rare signals. 

The rate of the newly observed mode B{B^ K-tt+) = (5.0 ± 0.75 ± 1.0) • 10"^ is 
in agreement with the latest theoretical expectation [29] which is lower than previous 
predictions [18, 30]. CDF measures for the first time in the B^ meson system the 
direct CP asymmetry Acp{B^s ^ K-7T+) = 0.39 ± 0.15 ± 0.08. This value favors 
a large CP violation in B^ mesons, on the other hand it is also compatible with 
zero. Ref. [16] suggests a robust test of Standard Model expectations versus new 
physics comparing the direct CP asymmetries in the S° — > K~7r'^ and B^ — > K'^n' 
decay modes. Using HFAG input [17], CDF measures r(B°^i^-7r+)-r(B0^i^+7r-) _ 

^ L J) r{BO^K-TT+)-r{B"^K+Tz-) 

0.84 ± 0.42 ± 0.15 (where F is the decay width) in agreement with the Standard 
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Model expectation of one. Assuming that the relationship above yields one and using 
as input the branching fraction — > K'tt'^) measured in this analysis, the world 
average for Acp{B^ K+n-) and the B{B^ ir+vr") [17], the expected value for 
Acp{B^s ~^ K^Ti^) ^ 0.37 is estimated in agreement with the CDF measurement. 

The rate of the mode K+R-) = (24.4 ± 1.4 ± 4.6) • 10"^ is in agreement 

with the latest theoretical expectation [31, 32] and with the previous CDF measure- 
ment [25]. An improved systematic uncertainty is expected for the final analysis of 
the same sample. The results for the meson are in agreement with world av- 
erage values [17]. The measurement Acp{B^ ^ K+tt') = -0.086 ± 0.023 ± 0.009 
is the world's second best measurement and the significance of the new world av- 
erage Ac^{B° K+n-) = -0.095 ± 0.013 moved from 6a to 7a. CDF updates 
the upper limits and quotes also the absolute branching fractions of the currently 
unobserved annihilation- type modes: B'^ K^K~ and B^ ti^ti~ . The rate 
B{B^ K+K-) = (0.39 ± 0.16 ± 0.12) ■ 10"^ has the same uncertainty as the cur- 
rent measurements [17], while the -B° tt+tt^ upper limit (already the world's best 
limit [25]) is improved by a factor of 1.3, approaching the expectations from recent 
calculations [19, 33] . CDF also reports the first observation of two new baryon charm- 
less modes pir^ and A° pK~ , and measures <B(A° p7r~)/B(A^ — > pK~) — 
0.66 ± 0.14 ± 0.08 in agreement with expectations from Ref. [34]. 

5 Summary 

We review recent result on heavy quark physics focusing on Run II measurements 
of B hadron spectroscopy and decay at the Tevatron. A wealth of new B physics 
measurements from CDF and DO has been available. These include the spectroscopy 
of excited B states {B**, B**) and the observation of the S5 baryon. The discus- 
sion of the decays of B hadrons and measurements of branching fractions focuses on 
charmless two-body decays of — > h'^h~ . We report several new 5° and A^ decay 
channels. 
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